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Long-duration (> 2 s) 7 -ray bursts that are believed to originate from the death of massive stars 
are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova¬ 
like emission down to very stringent limits and its physical origin is still debated. Here we report 
the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This 
red bump is inconsistent with even the weakest known supernova. However, it can arise from a 
Li-Paczyriski macronova — the radioactive decay of debris following a compact binary merger. If 
this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the 
death of a massive star and it was a site of a significant production of heavy r-process elements. The 
significant ejected mass favors a black hole-neutron star merger but a double neutron star merger 
cannot be ruled out. 


Long-duration (> 2 s) 7-ray bursts (GRBs) are believed to originate from Collapsars that involve death of massive 
stars and are expected to be accompanied by luminous supernovae (SNe). GRB 060614 was a nearby burst with a 
duration of 102 s at a redshift of 0.125 [H. While it is classified as a long burst according to its duration, extensive 
searches did not find any SNe-like emission down to limits hundreds of times fainter than SN 1998bw Q, the archetypal 
hypernova that accompanied long GRBs Moreover, the temporal lag and peak luminosity of GRB 060614 fell 

entirely within the short duration subclass and the properties of the host galaxy distinguish it from other long-duration 
GRB hosts. Thus, GRB 060614 did not fit into the standard picture in which long duration GRBs arise from the 
collapse of massive stars while short ones arise from compact binary mergers. It was nicknamed the “long-short burst” 
as its origin was unclear. Some speculated that it originated from compact binary merger and thus it is intrinsically a 
“short” GRB [l|, 0,0-0| • Others proposed that it was formed in a new type of a Collapsar which produces an energetic 
7—ray burst that is not accompanied by an SNe 

Two recent developments may shed a new light on the origin of this object. The first is the detection of a few 
very weak SNe (e.g. SN 2008ha i) with peak bolometric luminosities as low as L ~ 10^^ erg s The second is the 
detection of an infrared bump, again with a L ^ 10"*^ erg s“^, in the late afterglow of the short burst GRB 130603B [13, 
11|. This was interpreted as a Li-Paczyhski macronova (also called kilonova) |12 h 19I| —a near-infrared/optical transient 


powered by the radioactive decay of heavy elements synthesized in the ejecta of a compact binary merger. Motivated 
by these discoveries, we re-examined the afterglow data of this peculiar burst searching for a signal characteristic to 
one of these events. 

The X-ray and UV/optical afterglow data of GRB 060614, were extensively examined in the literature [13, [2l| and 
found to follow very well the fireball afterglow model up to t ~ 20 days [23|. The J-band has been disregarded because 
only upper limits ~ 19 — 20^^^ mag with a sizeable scatter are available at < > 2.7 day, and these are too bright to 
significantly constrain even supernovae as luminous as SN 1998bw [23|. Here we focus on the optical emission. We 
have re-analyzed all the late time (i.e., t > 1.7 day) Very Large Telescope (VLT) V, R and /—band archival data and 
the HST F606W and F814W archival data, including those reported in the literature [3,0| and several unpublished 
data points. Details on data reduction are given in the Appendix. Fig[T] depicts the most complete late-time optical 
light curves (see TabUin the Appendix; the late VLT upper limits are not shown in Figll} of this burst. 

The VLT V, R and /—band fluxes decrease with time as oc i-2.30±0.03 ^-ggg FiglU in which the VLT V/I band data 
have been calibrated to the F606W/F814W filters of HST with proper fc—corrections), consistent with that found 
earlier [3, [23[Ill. However, the first HST F814W data point is significantly above the same extrapolated power-law 
decline. The significance of the deviation is ~ 6 ct (see the estimate in the Appendix). No statistically-significant 
excess is present in both the F606W and the R bands. The F814W-band excess is made most forcibly by considering 
the color evolution of the transient, defined as the difference between the magnitudes in each filter, which evolves 
from V —/ ss 0.65 mag by the VLT (correspondingly for HST we have F606W — F814W ss 0.55 mag) at about t ~ 1.7 
day to F606W — F814W ~ 1.5 mag by HST at about 13.6 day after the trigger of the burst. With proper/minor 
extinction corrections, the optical to X-ray spectrum e nerg y distribution for GRB 060614 at the epoch of ~ 1.9 
day is nicely fitted by a single power-law /b oc ® [3, IS HH. In the standard external forward shock afterglow 
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FIG. 1: The afterglow emission, not corrected for the small amount of foreground and host extinction, of GRB 060614. Note 
that the VLT V/I band data have been calibrated to the HST F606W/F814W filters with proper fc—corrections (see the 
Appendix). The VLT data (the circles) are canonical hreball afterglow emission while the HST F814W detection (marked 
in the square) at t ~ 13.6 day is significantly in excess of the same extrapolated power-law decline (see the residual), which 
is at odds with the afterglow model. The F814W-band lightcurve of SN 2008ha expected at z = 0.125 is also presented for 
comparison. The dashed lines are Macronova model light curves generated from numerical simulation for the ejecta from 
a black hole—neutron star merger. 


model, the cooling frequency is expected to drop with time as Vc oc [i^. Thus, it cannot change the optical 

spectrum in the time interval of 1.9 — 13.6 day. Hence, the remarkable color change and the F814W-band excess of 
~ 1 mag suggest a new component. Like in GRB 130603B this component was observed at one epoch only. After 
the subtraction of the power-law decay component, the flux of the excess component decreased with time faster than 
^-3.2 Jqj. ^ ^ 23.6 days An unexpected optical re-brightening was also detected in GRB080503, another ‘long-short’ 
burst (^ . However, unlike the excess component identified here, that re-brightening was achromatic in optical to 
X-ray bands and therefore likely originated by a different process. 

Shortly after the discovery of GRB 060614 it was speculated that it is powered by an “unusual” core collapse of a 
massive star We turn now to explore whether the F814W-band excess can be powered by a weak supernova. 

Figl2]depicts the color F606W—F814W of the excess component (we take F606W—F814W~ 1.5 mag as a conservative 
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FIG. 2: The color (i.e., F606W-F814W) change of SN 2006aj, SN 2008ha, SN 2010bh and the “excess component” identihed 
in this work. The emission of these three supernovae, adopted from the literature [25I - [T7| has been shifted to 2 = 0.125, the 
redshift of GRB 060614, with corrections on the time, frequency and extinction. Note that the “excess component” is much 
redder than them. 


lower limit of the color of the “excess 
of SN 2006aj Q, SN 2010bh 


27| 


component due to the lack of simultaneous excess in F606W-band) with that 
and SN 2008ha, an extremely dim event [1^. The excess component has a 
much redder spectrum than the three supernovae. If the “excess component” was thermal it had a low effectively 
temperature < 3000 K to yield the very soft spectrum. Such unusually low effective temperature is also needed 
to account for the very rapid decline of the excess component. The expansion velocity can be estimated as u ^ 

“^(</13.6 day)“^. The implied ®®Ni mass is ^ 10“^ Mq if this 

i. : 


1.2 X 10^ km s-i (L/lO^i erg s-i)i/2(reff/3000 K) 


was a supernova-like event that peaked at ~ 13.6 day [9j. In this work we take a standard cosmology model with 
Hq = 71 km Mpc“^, Hm = 0.32 and Da = 0.68. Comparing with the extremely-faint SN 2008ha after proper 
corrections to z = 0.125, the peak F814W-band emission of the “excess component” is lower by ~ 1 mag and the 
decline is also much faster. Hence the “excess component” is remarkably different from SN 2008ha. 

The low luminosity as well as the low effective temperature of the transient emission are typical characteristics 
of a Macronova, a transient arising from the radioactive /3-decay of material ejected in a compact binary merger. 
The opacity of the macronova material is determined by the Lanthanides that are produced via r—process in the 
neutron-rich outflow. This opacity is very large (k ~ 10 cm^ g~^) resulting in a weak, late and red emission. The 
emerging flux is greatly diminished by line blanketing, with the radiation peaking in the near-infrared and being 
produced over a timescale of ~ 1 — 2 weeks [MIl Simple analytic estimates, using a radioactive /3-decay heating 
rate of lO^o erg s'^ g-^ [f/(l + z)l day]-i nii13. suggest that in order to explain the observed F814W-band 
excess, the required ejecta mass and expansion velocity are: Mej ^ 0.13 Mq (L/10'^^ erg s“^)(t/13.6 day)^'^ and 
V ^ O.lc (L/IO'^^ erg /2000 K)“^(t/13.6 day)“^, respectively. Note that the Macronova outflow is quite 

cold at such a late time [l7|, [l8| . The effective temperature is ~ 2000 K and the observer’s F814W-band is above 
the peak of the black body spectrum. The emitting radius and the corresponding expansion velocity are much larger 
than in a supernova at this stage. Scaled up numerical simulations of lighter ejecta from black hole—neutron star 
mergers suggest that Mej ~ 0.1 Mq and a velocity ^ 0.2c can account for the observed F814W-band excess. This 
numerical example is presented in Fig IT] in dashed lines. 

The implied ejecta mass is large compared with the mass ejection estimated numerically to take place in double 
neutron star mergers. However, it is within the possible range of dynamical ejecta of black hole-neutron star mergers 
with some extreme parameters (a large neutron star radius and a high black-hole spin aligned with the orbital angular 
momentum) 


An accretion disk wind may contribute some additional mass as well [15|, |33|, [3^ . However 


the radioactive heating due to hssion of the heavy r—process nuclei, which is quite uncertain and subdominant in 
current heating estimates [l^ . may play an important roleju the energy deposition. It may increase the energy 
deposition rate at around ten days by a significant factor 


35|. This may reduces the required ejecta mass to 


0.03 — 0.05 Mq. This range of the ejecta masses is well within the range of the dynamical ejecta of black-hole neutron 
star mergers and it is even compatible with some estimates of double neutron star mergers. 

We conclude that while a weak supernova cannot explain the observations, a high mass ejection macronova may. 
Like in GRB 130603B we must caution here that this interpretation is based on a single data point. However, if this 
interpretation is correct, it has far reaching implications. First, the presence of macronovae in both the canonical 
short burst GRB 130603B and in this “long-short” one, GRB 060614, suggests that the phenomenon is common and 
the prospects of detecting these transients are promising. A more conclusive detection based on more than a single 
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data point could be achieved in the future provided that denser HST observations are carried out. Moreover, as a 
black hole-neutron star merger is favored in explaining the large ejected mass this implies that such binary systems 
may exist and their mergers are also responsible for GRBs. It also suggests that the “long-short” burst was in fact 
“short” in nature, namely, it arose from a merger and not from a Collapsar. The fact that a merger generates a 100 
sec long burst is interesting and puzzling by itself. 

Clearly such events would contribute a significant fraction of the r—process material??. The actual contribution 
relative to the contribution of 130603B-like events is difficult to estimate as it is unclear which fraction of the 
Macronovae/kilonovae behave as each type. Because of beaming most mergers will not be observed as GRBs. 
However, they emit omnidirectional gravitational radiation that can be detected by the upcoming Advanced 
LIGO/VIRGO/KAGRA detectors. These near-infrared/optical macronovae could serve as promising electromagnetic 
counterparts of gravitational wave triggers in the upcoming Advanced LIGO/VIRGO/KAGRA era. 
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Appendix 


DATA REDUCTION 


We retrieved the public VLT imaging data of GRB 060614 from ESO Science Archive Facility 
(http://archive.eso.org). The raw data were reduced following standard procedures, including bias sub¬ 
traction, flat fielding, bad pixel removal, and combination. Observations made with the same instru¬ 
ment and filter at different epochs are compared to that of the last epoch. The software package ISIS 
(http://www2.iap.fr/users/alard/package.html) is used to subtract images and measure the GRB afterglow from 
the residual images. Photometric errors are estimated from the photon noise and the sky variance to Icr confidence 
level. The 3 ct of the background RMS of the residual images is taken as the limiting magnitude. Finally, standard 
stars observed on Jun 16, 2006 were used for the absolute calibration. The results are shown in TabUl We assumed 
that the afterglow is characterized by the same power-law spectrum with index j3 = 0.80 [ioj i during these observa¬ 
tions, with which we get the /c-corrections between the VLT V/I and HST F606W/F814W magnitudes, namely 0.12 
mag and 0.02 mag, respectively. Such corrections had been taken into account in FiglTJ 

HST archive data of GRB 060614 are available from the Mikulski Archive for Space Telescopes (MAST; 
http://archive.stsci.edu), including one observation with WFPC2 and four observations with ACS in F606W and 
F8I4W bands. The reduced data provided by MAST were used in our analysis. The last visit is taken as the 
reference and the other images of the same filter are subtracted in order to directly measure fluxes of the afterglow 
from the residual images. Empirical PSFs was built with bright stars in each image. Bright compact objects in the 
same field were used to align and relatively calibrate these images. WFPC2 image differs from ACS image in PSF. 
Before image subtraction, the WFPC2 and ACS images were matched to the same resolution by convolving each with 
the other’s PSF. The PSF-matched WFPC2 and ACS images were aligned and subtracted. Aperture photometry 
was carried out for the afterglow in the residual image. The aperture correction derived from the empirical PSF 
was applied to yield the total flux. The host galaxy was used to relatively calibrate the afterglow between images, 
and the ACS zeropoints were used for absolute calibration. If the signal of the afterglow is too faint to be a secure 
detection, an upper limit of 3 ct background RMS is adopted. The magnitudes of the host galaxy are measured in 
the last observation of all filters and can well be fitted by an Sc type galaxy template, see FiglU demonstrating the 
self-consistence of our results. Our results are summarized in TabU being well consistent with these given by other 
groups iSEl- 


THE DECLINE RATE OF THE VLT AFTERGLOW DATA AND THE SIGNIFICANCE OF THE EXCESS 


As found in previous studies, the late time optical/X-ray afterglow emission of GRB 060614 can be interpreted within 
the fireball forward shock model 2^ 211. Motivated by such a fact, we assume that the I, R, V lightcurves follow the 
same power-law decline. In our fit there are four free parameters, three are related to the initial flux/magnitude in 
these three bands and the last is the decline rate needed in further analysis. We fitted all the VLT data (combined 
I, R, V band together) during the first 15 days (after which there are just upper limits) to determine these four 
parameters as well as their errors. The best-fit decline is found to be oc consistent with that obtained 

in optical to X-ray bands in previous studies [1, [2^ [2lj. As a result of the propagation of uncertainties, the errors 
of the best-fit light curves are consequently inferred (the shadow regions in the residual plot of Fig.l represent the 
Icr errors of the best-ht light curves). Please note that in Fig.l the VLT V/I band emission have been calibrated to 
HST F606W/F814W filters with proper fc—corrections. The flux separation between the HST F814W-band data and 
the fitted curve at t ^ 13.6 day is Fexcess = 0.182 /rJy. The flux error of the F814W-band emission at < ~ 13.6 day 
is SFohs ~ 0.024 ^Jy. The flux error of the best fitted F814W-band lightcurve at f 13.6 days is SF^t ~ 0.012 ^Jy. 
The signihcance of the excess component is estimated by 7^ = Fe^cess/\/SF'^^^ + SF^^. ~ 6 . We therefore suggest that 
the excess component identified in this work is statistically significant at a confidence level of ~ 6tT. 
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TABLE I: Log of observations 


Time from GRB“ 
(days) 

Filter 

Exposure time Instrument 

(s) 

Vega Magnitude® 

1.72034 

V 

2 x120 

VLT+FORSl 

21.38+0.03 

2.83515 

V 

4x90 

VLT+FORSl 

22.88+0.07 

3.86077 

V 

2x120+4x180 VLT+FORSl 

23.64+0.09 

7.82790 

V 

3x180 

VLT+FORS2 

24.85+0.25 

23.79964 

V 

2 x120 

VLT+FORSl 

>24.8 

32.78644 

V 

5x120 

VLT+FORSl 

[22.96 + 0.04] 

108.57188 

V 

3x300 

VLT+FORS2 

[22.87 ± 0.04] 

1.72583 

R 

2 x120 

VLT+FORSl 

21.03+0.02 

1.86974 

R 

2 x120 

VLT+FORSl 

21.28+0.03 

2.84199 

R 

2 x120 

VLT+FORSl 

22.44+0.05 

3.86899 

R 

2x120+4x180 

VLT+FORSl 

23.14+0.04 

4.84365 

R 

2x180 

VLT+FORSl 

23.54+0.06 

6.74083 

R 

3x180 

VLT+FORSl 

24.28+0.09 

10.81441 

R 

2x300 

VLT+FORSl 

25.54+0.26 

14.77259 

R 

4x300+4x180 

VLT+FORSl 

26.35+0.32 

19.67818 

R 

6x240 

VLT+FORSl 

>26.3 

23.80494 

R 

2 x120 

VLT+FORSl 

>24.6 

32.79667 

R 

3x180 

VLT+FORSl 

>25.4 

44.73601 

R 

5x240+2x180 

VLT+FORSl 

>26.3 

64.70367 

R 

12x300 

VLT+FORSl 

[22.42 ± 0.03] 

1.73236 

I 

3x120 

VLT+FORSl 

20.73+0.04 

2.84826 

1 

3x120 

VLT+FORSl 

21.97+0.17 

3.85840 

I 

4x300 

VLT+FORSl 

22.48+0.10 

7.84052 

I 

3x120 

VLT+FORS2 

23.99+0.20 

23.81008 

I 

2 x120 

VLT+FORSl 

>23.9 

32.80572 
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Notes. 

a. Time since the burst trigger, which occurred at 2006 June 14, 12:43:48 UT. 
b. These values have not been corrected for the Galactic extinction of Ay = 0.07 mag. The magnitudes in square brackets are 
for the host galaxy. The definitions of the errors and upper limits are described in the text (note that 0.02 magnitude of 
uncertainty in absolutely calibrations has been added to the statistical errors). 
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FIG. 3: Following Gal-Yam et al. (2006)0, we fit the host galaxy magnitudes by an Sc type template. The redshift of the 
host galaxy 2 = 0.125 and the Galactic extinction Ay = 0.07 mag have been taken into account. 




















